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Abstract We present preliminary results on a study of the 2-850 SEDs of a sample 
of 30 FIRBACK galaxies selected at 170 /jm. These sources are representative 
of the brightest ~ 10% of the Cosmic Infrared Background. They are a mixture 
of mostly local (z < 0.3) starforming galaxies, and a tail of ULIGs that extend 
up to z~ 1, and are likely to be a similar population to faint SCUBA sources. 
We use archival Spitzer IRAC and MIPS data to extend the spectral coverage to 
the mid-IR regime, resulting in an unprecended (for this redshift range) census 
of their infrared SEDs. This allows us to study in far greater detail this im- 
portant population linking the near-IR stellar emission with PAH and thermal 
dust emission. We do this using a Markov Chain Monte Carlo method, which 
easily allows for the inclusion of " 6 free parameters, as well as an estimate of 
parameter uncertainties and correlations. 
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Introduction 

Understanding the full infrared spectral energy distributions (SEDs) of galax- 
ies is essential for a complete picture of star-formation in the Universe. About 
half the power ever generated by stars is reprocessed by dust into the infrared as 
inferred from the Cosmic Infrared Background (e.g. Puget et al. 1996; Dwek et 
al. 1998; Hauser et al. 1998). This emission is increasingly important at higher 
redshifts where the star-formation density of the Universe is larger than today, 
coupled with our observational bias toward detecting more active galaxies at 
high-z. The variation in SED shapes is a key uncertainty when comparing pop- 
ulations selected at different wavelengths or testing galaxy evolution models. 

We define the infrared SED as covering the range 3-1000 /im, which roughly 
spans the wavelengths between fully stellar and fully synchrotron-dominated 
emission. Traditionally, studying this entire range at once has been difficult, 
since mid-IR (< 60 mJy) observations could not reach much beyond the local 
Universe (except, to some extend, for the deepest ISOCAM 15 fim observa- 
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tions), while sub-mm observations also only exist for very local, IR-bright, 
galaxies (e.g. Dunne et al. 2000), or else SCUBA-selected sources which span 
z ~ 1-3 (Chapman et al. 2003; Pope et al. 2005). The latter typically have only 
one or two detections outside the single SCUBA (850 /im) one, which makes 
their interpretation particularly dependent on the SED model assumed (Blain 
et al. 2002). Due to these past observational limitations, we still do not know 
(beyond some generalized trends) the full range of galaxy SED shapes, how 
exactly they are related to the underlying physical conditions in the galaxy, 
and therefore how they may vary across cosmic time as the galaxies evolve. 

With the advent of the Spitzer Space Telescope (Werner et al. 2004) we can 
for the first time observe the mid-IR (3 - 24 fim) properties of large numbers of 
sources over a cosmologically significant range in z (e.g. Lonsdale et al. 2003). 
The obvious next step toward characterizing the full infrared SEDs of galax- 
ies is therefore to link Spitzer observations with longer wavelength samples, 
especially including sub-mm observations. 

The sample of 30 galaxies discussed here is a sub-sample from the 170 /im 
FIRBACK (Far-IR BACKground) ELAIS-N1 catalog (Dole et al. 2001). They 
were selected on the basis of an existing radio detection which we followed-up 
with both deep near-IR and sub-mm observations (Scott et al. 2000; Sajina 
et al. 2003). The sample selection appears to be unbiased with respect to 
the FIRBACK population as a whole. Our previous studies suggest that the 
sample consists primarily of z < 0.3 ordinary spiral-like galaxies rich in cold 
(T < 40 K) dust, with ~l/6 of the sample consisting of Ultraluminous Infrared 
Galaxies (ULIGs) at z~ 0.5-1. Spectroscopic follow-up of this and related 
sub-samples confirm both the local sources (Patris et al. 2003; Dennefeld et 
al. 2004), and the higher- z ones (Chapman et al. 2002). Thus this sample rep- 
resents a bridge population between the local Universe and SCUBA blank-sky 
sources (e.g. Pope et al. 2005). 

Here we use archival Spitzer observations of the ELAIS-N1 field (see Fig 1) 
in order to extend the known SEDs of the above sample into the mid-IR wave- 
length range. The IRAC and MIPS plus ISO photometry, together with ground- 
based near-IR and sub-mm data, gives us more than 10 constraints for every 
galaxy. We fit these SEDs with a phenomenological model motivated by differ- 
ent physical origins for the emission. We investigate statistically robust ways 
to fit the SEDs based on Markov Chain Monte Carlo (MCMC). These allow 
us to go beyond merely finding the best-fit SED model, toward studying the 
spread in multi-dimensional, inter-correlated SED models allowed by the data. 
This not only leads to reliable errors on the derived parameters (e.g. dust tem- 
perature, PAH emission strength), but also has obvious implications for the 
interpretation of data from current and future infrared surveys. 
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Figure 1. The brightest 12 of our sources as seen by Spitzer, where red, blue, and green are 
respectively the 8 fim, 3.6 fim, and 5.8 fim IRAC channels. The box size is 2' (roughly enclosing 
the ISO 170 pm beam). 

SED model & MCMC fitting 

We model the SED as a sum of stellar emission, PAH emission, power-law 
emission, and thermal grey-body emission. This is similar to the model used 
in Sajina, Lacy, and Scott (2005) except for the addition of the thermal com- 
ponent to account for the far-IR and sub-mm emission. The other modification 
is in the power-law component which now mostly accounts for the 12-100 /im 
emission and is cut-off at long-wavelengths in order not to contribute to the 
far-IR/sub-mm. The stellar emission is accounted for by a lOGyr-old, so- 
lar metallicity, Salpeter IMF, single stellar population (SSP) spectral template 
generated with the PEGASE2.0 spectral synthesis code (Fioc et al. 1997). This 
amounts to a six parameter model, although a seventh parameter (ry) is needed 
to fit the ULIG population where optical depths of ry ~ 10 are necessary. To 
find the best-fit model we use the genetic algorithm pikaia (Charbonneau et 
al. 1999). Once the loci of interest in the multi-dimensional parameter space 
have been found, we investigate the probability of the parameters via a Markov 
Chain Monte Carlo (MCMC) (see Sajina, Lacy, and Scott 2005 and references 
therein). Fig. 2 shows some examples of the model fits and their MCMC scat- 
ter. 
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Figure 2. Phenomenological modelling of the 2-850 fim SED of starforming galaxies. On 
the left, we show a representative sub-sample of our galaxies. The dotted line shows the cold 
dust component alone. The 8 pirn generally appears lower than the curve, since we account for 
the filter profiles. On the right, we show an example of the spread in allowable SED shapes 
from the MCMC fitting. The inset shows the marginalized, posterior probability distribution for 
the dust temperature. 



Discussion 

Here we have presented a new SED model, which aims to be a flexible, 
template-free approach for investigating the spectral properties of starforming 
galaxies. The Markov Chain Monte Carlo fitting allows us to sample the pos- 
terior probability distribution, providing reliable scatter in the derived parame- 
ters as well as revealing any parameter degeneracies. The spectral coverage of 
the sources shown here is comparable to that expected for upcoming surveys 
investigating the connections between the mid-IR and far-IR/sub-mm popula- 
tions (e.g. studies in the GOODS-North and SHADES fields). The advantage 
of this present sample, however, is that being largely local (z < 0.3), it displays 
a greater range in intrinsic properties (e.g. SFR), and we already know the red- 
shifts for most of the sources. This also allows us to test various photometric 
redshift estimators. Detailed analysis of the results of this fitting procedure, 
and implications for the FIRBACK population is work in progress (Sajina et 
al. 2005b). 
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